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We present first-principles calculations of the sticking coefficient of O 2 at Pd(lOO) to assess the effect of 
phononic energy dissipation on this kinetic parameter. For this we augment dynamical simulations on six¬ 
dimensional potential energy surfaces (PESs) representing the molecular degrees of freedom with various 
effective accounts of surface mobility. In comparison to the prevalent frozen-surface approach energy dissi¬ 
pation is found to qualitatively affect the calculated sticking curves. At the level of a generalized Langevin 
oscillator (GLO) model we achieve good agreement with experimental data. The agreement is similarly 
reached for PESs based on two different semi-local density-functional theory functionals. This robustness of 
the simulated sticking curve does not extend to the underlying adsorption mechanism, which is predominantly 
directly dissociative for one functional or molecularly trapped for the other. Completely different adsorption 
mechanisms therewith lead to rather similar sticking curves that agree equally well with the experimental 
data. This highlights the danger of the prevalent practice to extract corresponding mechanistic details from 
simple fingerprints of measured sticking data for such exothermic surface reactions. 

PACS numbers: 34.35.-ka,68.49.Df,82.20.Kh 


I. INTRODUCTION 

Accurate knowledge of the sticking coefficient of 
molecules at transition metal surfaces is generally valu¬ 
able in view of the relevance of this kinetic parameter in 
a wide range of applications such as heterogeneous catal¬ 
ysis. On more reactive surfaces (often desirable) dissoci¬ 
ation can already happen alongside the adsorption pro¬ 
cess. Sticking data may then further provide informa¬ 
tion on the underlying adsorption mechanism. On the 
other hand, a sticking coefficient is an averaged kinetic 
quantity, comprising contributions from initially varying 
translational and rotational molecular impingement. To 
what extent sticking data really reflects mechanistic de¬ 
tails is thus an open question^. Prom a modeling perspec¬ 
tive, there is vice versa the question of how well these un¬ 
derlying details need to be described to arrive at reliable 
sticking coefficients that can in turn e.g. be employed in 
microkinetic models^. 

Aiming to answer these questions, intense fundamental 
research has been conducted to obtain an atomic-scale 
understanding of the reaction pathways and dynamics 
governing the adsorption process. Supersonic molecu¬ 
lar beam techniques at low-index single crystal surfaces 
have made an unequivocal contribution in this respect, 
providing sticking coefficient data for a wide range of 
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well-defined initial conditions of incidence energy and an¬ 
gle, as well as resolved rotational or vibrational state^. 
In terms of dissociative adsorption corresponding data 
has been traditionally analyzed to distinguish regimes 
of direct or indirect adsorption^. The former essentially 
suggests immediate dissociation upon impact with the 
surface. The latter alludes to the notion of equilibration 
with the surface in a temporary intermediate (molecular) 
state which serves as a precursor to dissociation. Detailed 
studies of prototypical diatomics have especially served 
to establish general trends for sticking curves, i.e. certain 
signatures are seen to indicate one or the other adsorption 
mechanism®. A dependence of sticking data on substrate 
temperature is for instance believed to reflect precursor- 
mediated adsorption, considering the kinetic competition 
between dissociation and desorption from such a precur¬ 
sor state. 

Such trends have met an overall intriguing consistency 
over a range of studied systems. Notwithstanding, one 
needs to recognize that a bulk of corresponding work 
has been done for (dissociative) H 2 adsorption, where a 
clearcut distinction might be facilitated by the relatively 
“smooth” molecule-surface interaction. Upon presence 
of TT-orbital involving molecular bonds, the correspond¬ 
ing potential energy surface (PES) becomes significantly 
more intricate, which could prohibit such a direct ex¬ 
traction of mechanistic insight from measured sticking 
data. A characteristic example for this is the adsorp¬ 
tion of N 2 on W(llO), where measured sticking curves 
eluded a simple interpretation on the basis of the afore- 



2 


mentioned experimental trends®. In fact, seemingly con¬ 
tradictory evidence from molecular beam measurements 
could only be reconciled through explicit dynamical sim¬ 
ulations, which revealed a complex interplay between a 
direct and indirect adsorption mechanism depending on 
initial conditions®. The contribution of these processes 
was shown to arise from characteristics of the gas-surface 
interaction far from the surface and could only be cap¬ 
tured within an accurate high-dimensional description of 
the PES involving all molecular degrees of freedom. 

Similar complexities can be expected for the adsorp¬ 
tion of O 2 , a key process in oxidation catalysis. In fact, 
the typically large exothermocity of this process on tran¬ 
sition metal surfaces adds yet another facet, namely the 
question of energy dissipation. Considering its influence 
e.g. on the aforementioned kinetic competition of dissoci¬ 
ation or desorption from a precursor state, the question 
is how much this affects measured sticking curves and 
the established trends in terms of fingerprints. From a 
modeling perspective it again adds the question of how 
well the dissipative dynamics needs to be accounted for. 
The latter is particularly relevant, as a fully quantita¬ 
tive account of phononic energy dissipation in explicit 
ab initio molecular dynamics (MD) simulations is still 
highly demanding, when aiming for sufficient statistical 
averaging and the full computation of energy-dependent 
sticking curves^’®. Highly appealing are therefore more 
effective treatments of surface mobility®^^®. The poten- 
tiaf sensitivity of O 2 sticking curves to the details of en¬ 
ergy dissipation then makes a comparison to high-quality 
experimental data particularly valuable to gauge the ac¬ 
curacy of such effective approaches or vice versa the level 
of detail required to account for. 

With this motivation we focus in the present work on 
the adsorption of O 2 on clean Pd(lOO). Recent molecu¬ 
lar beam experiments found the initial sticking probabil¬ 
ity So{E[,Ts) for this system to be independent of sub¬ 
strate temperature Tg and only weakly increasing with 
incident kinetic energy In terms of the classical 

trends, this suggested an interpretation in form of a pre¬ 
dominantly direct dissociation mechanism with, at most, 
some reaction paths that include a modest activation 
barrier. At low Ei and Tg, however, the contribution 
from a partly equilibrated molecular-precursor was con¬ 
jectured in order to rationalize the independence on sur¬ 
face coverage and deviation from normal energy scaling 
under those conditions'^. We scrutinize this interpreta¬ 
tion through dynamical simulations on a first-principles 
six-dimensional (6D) PES that accounts for all molecu¬ 
lar degrees of freedom and which we suitably augment 
with effective treatments of surface mobility. Accounting 
in some respects for the latter is found to substantially 
change the calculated sticking curve So{Ei,Ts), i.e. the 
latter is indeed sensitive to energy dissipation. Intrigu- 
ingly this holds for both direct and indirect adsorption 
mechanisms, either of which we obtain as dominant when 
basing the simulations on PESs obtained with two dif¬ 
ferent density-functional theory (DFT) functionals. The 


uncertainties introduced by current semi-local DFT func¬ 
tionals thus prohibit a clear identification of the dom¬ 
inant adsorption mechanism. They also do not allow 
to fully disentangle whether the approximate treatment 
of substrate mobility or the deficiencies in the underly¬ 
ing DFT energetics are the primary reason for remaining 
small differences to the experimental data. All these in¬ 
tricacies nevertheless point already at this stage at the 
limitations of trying to directly deduce insight on the 
character of the adsorption process from measured stick¬ 
ing curves alone. 


II. METHODS 

The O2-Pd(100) interaction energetics are obtained 
by spin-polarized DFT calculations, using either the ex¬ 
change correlation functional due to Perdew, Burke and 
Ernzerhof (PBE)^®d6 qj. Hammer, Hansen and 

Nprskov (RPBE)^^ to approximately assess the uncer¬ 
tainties introduced by prevalent semi-local functionals 
(see below). Electronic states are described with a plane 
wave basis set using a cut-off energy of 400 eV as im¬ 
plemented in the CASTEP code^®, together with ultrasoft 
pseudopotentials (USPPs)^® as bundled in the Materi¬ 
als Studio 6.0 database. These USPPs have been ob¬ 
tained with Vanderbilt’s original generator^® using the 
PBE functional. By comparing to a few converged all¬ 
electron calculations based on the FHI-aims code^®, we 
have verified that the pseudopotential-induced error^^’^^ 
in the O2-Pd(100) interaction energies is not larger than 
100 meV^®. Within a periodic supercell model, the sur¬ 
face is represented by five-layer slabs which are separated 
by a vacuum distance of 15 A and which are (3x3) multi¬ 
ples of the primitive surface unit cell of Pd(lOO). Calcu¬ 
lations are performed using a (4x4x1) Monkhorst-Pack 
grid^^ for k-point sampling. 

More than 6,000 DFT energies are calculated for var¬ 
ious high- and low-symmetry configurations of the oxy¬ 
gen molecule above the frozen Pd(lOO) surface. These 
provide the basis for constructing a continuous repre¬ 
sentation of the adiabatic PES within all six molecular 
degrees of freedom (V6 d) for each of the two DFT func¬ 
tionals based on symmetry-adapted neural networks as 
pioneered by Behler and Reuter^® for fcc(lll) surfaces. 
Details about the present implementation can be found 
in Ref.^®’^®’^^, including an adaption for fcc(lOO) surfaces 
in particular. We note that the resulting DFT-PBE PES 
has already been employed previously in Ref. 8, 23, and 
28. Exactly the same recipe is followed here to obtain 
the continuous PES within the RPBE functional. See 
supplemental material for more details about the qual¬ 
ity of both NN PES representations. Minima and bar¬ 
rier searches on these PESs have been conducted with a 
stochastic sampling method described in earlier work^^ 
and with the Nudged Elastic Band (NEB) method as im¬ 
plemented within the Atomic Simulation Environment 
(ASE)^®, respectively. 
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Classical MD simulations are performed on the result¬ 
ing numerically efficient continuous PES representations. 
The effect of including the initial zero point energy of 
the O 2 molecule within a quasi-classical treatment was 
additionally investigated but found to have a negligi¬ 
ble effect on the resulting dynamics. Focusing only on 
the molecular degrees of freedom, such simulations do 
not allow for phononic energy dissipation and will hence¬ 
forth be denoted with frozen surface (FS). A first account 
of surface mobility can be incorporated on the level of 
the 3D surface oscillator (SO) model^. Here the sur¬ 
face is mimicked by an oscillator which is assigned the 
mass (toso) of a single Pd atom and is permitted to 
move as a whole in all three directions within a har¬ 
monic potential. The associated (3 x 3) frequency matrix 
ihso is assumed to be diagonal with values correspond¬ 
ing to a well localized surface mode of the Pd surface®’^^: 

= ^wsOyy = 16meV and huiso^^ = llmeV. The 
02-phonon coupling is then described by a 3D-space rigid 
shift Rso = (Aso,bsO) ^so) of VeD, and the MD equa¬ 
tions of motion in this approximation are given by 


9^R-a,b 


S^Rso 


-bF6d(Ra - Rso; Rb — Rso) 

mA,B 

(la) 

-VrsoF6d(Ra — Rso; Rb — Rso) 

mso 

- WgQ ■ Rso , (lb) 


where mA/e and Ra/b are the masses and Cartesian co¬ 
ordinates of the two individual oxygen atoms A and B. 


The effect of a bulk thermal bath is approximately in¬ 
cluded within the generalized Langevin oscillator (GLO) 
approach^°“^^. Here, the SO is coupled to a further 3D 
so-called ghost oscillator of equal mass, moLO = ''tt-SOi 
and frequency matrix wglo = wso- The same fre¬ 
quencies are also used to describe the SO-GLO coupling 
through a (3 X 3) coupling matrix Aso-glo = ^SO- 
The ghost oscillator is subject to frictional and random 
forces in order to account for energy dissipation and ther¬ 
mal fluctuations, respectively. As originally proposed by 
Adelman and Doll^*’, the former are described through an 
isotropic and diagonal damping matrix 7glo = 7glo 1 = 
ttwd/S 1 , where ujb is the Pd bulk Debye frequency^®’®®. 
Finally, the random force is a Gaussian white noise source 
W with a variance of (2fcBrs7GLo/”^GLO^0^^^i where 
/cb is the Boltzmann constant and At is the MD time in¬ 
tegration step. The resulting equations of motion within 


the GLO model are thereby as follows: 
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We note that this GLO implementation and specific 
choice of parameters follows that of earlier works regard¬ 
ing H 2 adsorbing on or scattered from the Pd(lll) and 
Pd(llO) surfaces®^^®®. We nevertheless systematically 
tested the dependence on the specihc parameter values by 
varying the oscillation frequencies entering urso, <^GhO, 
and Aso-glo by one order of magnitude, by varying the 
damping coefficient 7glo by two orders of magnitude, as 
well as by doubling the mass mso- This had little ef¬ 
fect on the simulation results as will be further specified 
below. 

The initial sticking coefficient at normal O 2 incidence 
was determined from classical MD trajectories with the 
O 2 molecule initially with its center of mass at a distance 
Z = 9 A from the surface, where the PES value for oxy¬ 
gen at its equilibrium bond length d^q = 1.24 A is zero 
in both the DFT-PBE and DFT-RPBE PES representa¬ 
tions. The initial molecular orientation and lateral cen¬ 
ter of mass position were sampled using a conventional 
Monte Carlo procedure. All statistical quantities are ob¬ 
tained by averaging over 5,000 trajectories for each value 
of incidence energy Ei and substrate temperature Tg. In¬ 
dividual trajectories were integrated up to 10 ps in order 
to reach the following classification in terms of molecular 
adsorption, dissociative adsorption or reflection: A tra¬ 
jectory was classified as dissociative whenever the O 2 in- 
ternuclear distance d reached twice its equilibrium value 
{d > 2deq) and is further increasing at this time (d > 0), 
while reflection was concluded when the molecular cen¬ 
ter reached its initial starting distance above the surface 
with a positive Z-velocity. All trajectories where neither 
dissociation nor reflection occurred up to the 10 ps inte¬ 
gration time were classified as trapped. 


III. RESULTS AND DISCUSSION 

A. Experimental sticking versus dynamics within the 
frozen surface approximation 

Figure 1 shows the initial sticking probability 
S'o(Fi,Ts) at normal incidence as recently measured by 
molecular beam experiments^^. Sticking is generally high 
(> 70%) within the entire range of investigated inci¬ 
dence energies and largely identical at the two substrate 
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Incidence energy, Ei (eV) 


FIG. 1. (Color online) Initial sticking probability So{Ei, Tb) of 
O 2 on Pd(lOO) as a function of incidence energy Ei and at nor¬ 
mal angle of incidence. Experimental data for two substrate 
temperatures, Tb = 100 K (solid blue line) and 400 K (solid red 
line), is reproduced from Ref. 14. Theoretical sticking proba¬ 
bilities are calculated in the temperature-independent frozen- 
surface (FS) approximation, using either DFT-PBE (dashed 
black line) or DFT-PRBE (dotted black line) energetics. 


temperatures investigated, 100 K and 400 K. The overall 
weak increase with incidence energy and independence 
of substrate temperature was interpreted to reflect pre¬ 
dominantly direct dissociative adsorption^^. From addi¬ 
tional data taken at finite coverages and oblique angles of 
incidence, some contribution from a partly equilibrated 
molecular-precursor was nevertheless speculated at low 
(EuTb). 

Even without explicit dynamical simulations the over¬ 
all high reactivity can already be gleaned directly from 
the calculated O2-Pd(100) interaction potential. From 
a global search on the continuous 6D DFT-PBE PES 
representation^^, we identify as most stable molecular ad¬ 
sorption state a configuration in which the O 2 molecule 
centers side-on above the Pd(lOO) hollow site with its 
molecular axis oriented along the [001] direction, i.e. in 
the direction of the neighboring hollow sites. This molec¬ 
ular well can be reached along a barrierless entrance 
channel as visible in the 2D (elbow) cut through the 
PES shown in Fig. 2. The calculated minimum energy 
path to dissociation also included in Fig. 2 exhibits a 
dissociation barrier out of this well of Ediss = 200 meV. 
This is a small value compared to the actual depth of 
the molecular well of Eb = — 1.45eV with respect to the 
gas phase. Even if a fraction of this energy gained upon 
adsorption is quickly transferred to internal molecular 
degrees of freedom, such a small barrier should thus be 
readily surmounted. Based on highly consistent DFT- 
PBE energetics (Eb = —1.52eV, E^iss = 120meV), this 
was the conclusion taken by Liu and Evans, who cor¬ 
respondingly classified the dissociative adsorption pro¬ 
cess as essentially non-activated^^. The overall high re¬ 
activity measured in the experiment seems to support 



FIG. 2. (Color online) Contour plot of a 2D cut through the 
6D O2/Pd(100) DFT-PBE PES depicting the existence of a 
pronounced molecular chemisorption well. Shown is the en¬ 
ergy profile as a function of the vertical distance of the O 2 
center of mass Z and the O 2 bond length d for a O 2 molecule 
in a side-on configuration with the molecular axis oriented 
along the [001] direction. Additionally shown by black mark¬ 
ers is the minimum energy path towards dissociation as ob¬ 
tained from a NEB calculation. Black contour lines indicate 
250 meV energy increments. Inset: Lateral O 2 center of mass 
positions (X, Y) over the surface unit-cell for the ultimately 
dissociating trajectories as molecules with Ei — 25 meV first 
reach a surface distance of ^ = 1.5 A (see text). 


this assignment, as well as the general notion to assume 
an approximately unity sticking coefficient in coarse¬ 
grained microkinetic simulations whenever static total 
energy calculations identify a barrierless entrance chan¬ 
nel to adsorption^’^®^^^. 

Aiming to scrutinize such mechanistic assignments 
made on the basis of the static PES alone we proceed 
with the classical trajectory calculations. Evaluating the 
fraction of reflected and adsorbing trajectories yields the 
theoretical sticking coefficient at a given incidence energy. 
Within the frozen-surface approximation, i.e. restriction 
to the 6D PES, the resulting sticking curve So{Ei) is 
necessarily independent of substrate temperature and is 
compared to the experimental sticking curves in Fig. 1. 
Confirming the expectations from the attractive DFT- 
PBE PES the theoretical sticking is rather high at low 
Ei and therewith in the same ballpark as the experimen¬ 
tal data. For these low incidence energies there is thus 
no qualitative disagreement with experiment at the level 
of semi-local DFT as e.g. in case of the enigmatic O 2 at 
Al(III) system^®’^®. This is not too surprising consider¬ 
ing that a major break-down of the electronic adiabatic- 
ity of the adsorption process had been invoked as one 
possible reason behind the 02 /Al(III) discrepancy®®’^®. 
This unlikely applies to Pd(IOO), which exhibits a very 
high electronic density-of-states at the Fermi level and 























5 


for which only small energy losses to electron-hole pairs 
during O 2 adsorption have been calculated^®. 

In terms of overall (rough) magnitude the favorable 
agreement of experimental and theoretical sticking ex¬ 
tends over the entire range of incidence energies shown 
in Fig. 1. Nevertheless, there is a disturbing difference in 
trend with Ei. In contrast to the measured slight rise of 
sticking with Ei, the calculated FS sticking curve shows 
a continuous decline. This decline is in fact even more 
pronounced when calculating the sticking curve on the 
basis of the DFT-RPBE PES, cf. Fig. 1. Repeating 
the above static analysis also for this PES we hnd its 
topology to be qualitatively very similar to that of the 
DFT-PBE PES. Important quantitative differences are 
therefore nicely summarized in terms of the two dynam¬ 
ically relevant quantities, Eb and E^iss- In line with the 
general construction idea of the RPBE functionab^, the 
molecular well is much more shallow in the DFT-RPBE 
PES (Eb = —0.85eV), while reflecting the concomitant 
weaker bond activation the dissociation barrier is with 
Ediss = 400 meV about twice as high as at the DFT- 
PBE level. Even though the ratio of these two quantities 
is thus much less favorable, on purely energetic grounds 
dissociative adsorption would nevertheless still be classi¬ 
fied as non-activated even at the DFT-RPBE level, i.e. 
the molecules gain much more energy upon adsorption 
than is needed to overcome the dissociation barrier. 

The decline of the sticking curve with E, obtained with 
both functionals is thus a purely dynamical effect. Aim¬ 
ing to extract its origin we analyze the trajectories in 
more detail. First of all, this trajectory analysis largely 
confirms the dynamical relevance of the minimum en¬ 
ergy path depicted in Fig. 2. Independent of the in¬ 
cidence energy and for both functionals (as well as all 
later surface mobility treatments), essentially all trajec¬ 
tories leading to adsorption show the molecules accumu¬ 
lating first around the hollow molecular chemisorption 
well. This preference for adsorption above hollow is ex- 
emplihed in the inset of Fig. 2 which analyzes the lateral 
O 2 center of mass positions over the surface unit-cell for 
the ultimately dissociating dynamical trajectories as the 
molecules first reach a surface distance of Z = 1.5 A. 
The data shown corresponds to Ei = 25meV and the 
PBE functional, with equivalent findings obtained at all 
other incidence energies and for DFT-RPBE. 

After arrival at the molecular chemisorption state, the 
fate of the molecules is almost instantaneously decided 
at low Ei at the DFT-PBE level. This is demonstrated 
by the time distribution shown in Fig. 3a, which indi¬ 
cates after which time a trajectory has been classified as 
dissociated, reflected or trapped. Recall that trapping 
is assigned to all molecules that have neither dissociated 
nor reflected after 10 ps simulation time, which is why the 
corresponding fraction is shown at this time bin in Fig. 
3. As apparent from Fig. 3a, essentially all trajectories 
dissociate in a rather direct fashion, i.e. within several 
picoseconds and thus after a minimal number of surface 
rebounds (typically < 5). This is completely different at 
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FIG. 3. (Color online) Time distribution characterizing 
events of reflection (blue), direct dissociation (shaded red) 
and molecular trapping (gray) as classified within the per¬ 
formed MD simulations. Trapping is assigned if neither dis¬ 
sociation nor reflection has occurred after 10 ps simulation 
time, which is why the corresponding fraction is shown at 
this time bin. Left and right panels refer to trajectories cal¬ 
culated on the basis of the DFT-PBE and DFT-RPBE PES 
descriptions, respectively, while the employed level of theory 
and initial conditions {Ei and Ts) are as reported in the cor¬ 
responding panel labels. We note that summing up the red 
and gray bars yields the initial sticking probabilities So of the 
respective functional at the particular incidence energies Ei 
and surface temperatures Ts shown in Figs 1 and 4. 


the DFT-RPBE level, where almost all trajectories in¬ 
stead end up being trapped, cf. Fig. 3b. The different 
reactivity of the two DFT PESs {E^ and Ediss) thus does 
have a significant effect in terms of entirely changing the 
dominant adsorption mechanism from direct (DFT-PBE) 
to indirect (DFT-RPBE). At low incidence energies this 
does not show up in the sticking coefficient though, as 
already a small amount of energy transfer to internal vi¬ 
brational, rotational or lateral translational degrees of 
freedom is sufficient to prevent immediate desorption, 
and since molecules correspondingly trapped in the sur- 
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face potential are being counted as contributing towards 
sticking. We note that in the absence of any other dis¬ 
sipation mechanism within the FS approximation, such 
trapped molecules would in fact eventually be able to des¬ 
orb. Notwithstanding, the time scale for a correspond¬ 
ing transfer of energy back into the perpendicular trans¬ 
lational motion would be rather long. In reality other 
dissipation channels would have set in during this time, 
which is why it is reasonable to count these trajectories 
as contributing to the sticking coefficient. Also the rather 
arbitrarily set maximum integration time of 10 ps should 
in this respect not matter, as we have verified by applying 
the trapping classification already after 8 ps. 

While the PES differences thus do not show up at low 
Ei, they do increasingly at larger incidence energies. Cor¬ 
responding fast molecules are less efficiently steered to 
direct dissociation. Neither is the energy transfer to in¬ 
ternal degrees of freedom efficient enough anymore to 
quickly remove the large excess kinetic energy and trap 
them in the surface potential. As shown in Figs. 3c and 
d for the largest Ei = 425 meV we thus find for both 
PES descriptions a large number of molecules that are 
reflected very quickly from the surface. Without efficient 
trapping, only the decreasing fraction of molecules that 
already starts from configurations favorable for a more or 
less direct dissociation mechanism can contribute to the 
decreasing sticking seen in Fig. 1. Such a mechanism is 
generally disfavored by a smaller acceleration into a more 
shallow adsorption well and a higher barrier to surmount, 
which is why the concomitant decline of So(Ei) with Ei 
is also much more pronounced at the DFT-RPBE level. 

From a bulk of work on O 2 adsorption on late tran¬ 
sition metals^°, the PBE and RPBE functionals can be 
seen as popular representatives for opposite ends within 
the range of current gradient-corrected functionals, with 
the prior likely more on the overbinding side and the 
latter possibly slightly underbinding. Within the under¬ 
standing of the just presented dynamical analysis con¬ 
comitant uncertainties in the energetic description of the 
molecular chemisorption well {Ey, and E^iss) at the semi¬ 
local DFT level are nevertheless unlikely to cause the 
wrong trend of the sticking curve with Ei as compared 
to the experimental data. Rather than the overall attrac¬ 
tiveness of the PES, the major reason for the decline of 
the calculated sticking curves lies in an insufficient abil¬ 
ity to trap the molecules in the chemisorption well. This 
thus points more to a weakness of the FS approximation, 
in which such a trapping can only result from energy 
transfer into molecular rovibrational degrees of freedom. 
In reality, this can also come from inelastic collisions with 
the surface atoms, which suggests the necessity to include 
some degree of surface mobility into the modeling. 


B. Surface mobility on the level of SO and GLO models 

We introduce a first account of energy exchange with 
the lattice through the SO model and show the derived 
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FIG. 4. (Color online) Comparison of the initial sticking prob¬ 
ability So{Ei,Ta) of O 2 on Pd(lOO) as calculated within the 
surface oscillator model (SO) for normal angle of incidence 
and varying incidence energy Ei. Data are shown for substrate 
temperatures of E = 100 K (blue circles) and 400 K (red tri¬ 
angles), using either DFT-PBE (dashed lines) or DFT-PRBE 
(dotted lines) energetics. 


initial O 2 /Pd( 100 ) sticking probabilities for substrate 
temperatures of Tg = 100 K and 400 K in Fig. 4. On 
the basis of the DFT-PBE energetics, and compared to 
the equivalent description within the FS approximation, 
we indeed observe the qualitative change expected from 
our preceding dynamical analysis, namely a considerable 
increase in high-energy sticking (up to as much as 25% for 
Ei = 425 meV). This almost completely lifts the decline 
with Ei seen at the FS level. The resulting essentially 
constant So{Ei,Ts) is furthermore only weakly depen¬ 
dent on Ts and thus in seemingly improved agreement 
with experiment. 

The mechanistic analysis of the PBE-based adsorp¬ 
tive trajectories reveals the relevance of the same reac¬ 
tion pathway above hollow as found in the FS dynam¬ 
ics and depicted in Fig. 2. Additionally, the short reac¬ 
tion/reflection times analyzed in Fig. 3e once again sug¬ 
gest that the fate of impinging molecules is almost instan¬ 
taneously determined and direct dissociation takes place 
essentially upon first impact with the surface. The dif¬ 
ferences in reactivity compared to the FS case do there¬ 
fore not originate from the activation of different reaction 
pathways or adsorption mechanisms. We instead ratio¬ 
nalize the increased SO sticking at high Ei with the im¬ 
proved efficiency of direct dissociation above hollow due 
to an enhanced steering of impinging fast molecules in the 
now “elastic” surface potential. The energy transfer to 
the surface possible within the SO model simply damp¬ 
ens the perpendicular translational motion of highly ac¬ 
celerated O 2 already upon first impact. This increases 
the probability that they will be effectively captured by 
the attractive potential and subsequently led to dissoci¬ 
ation. This picture arises from the analysis of the sig¬ 
nificant amount of trajectories that start from the exact 
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same initial conditions and yet are reflected within the FS 
but dissociated within the SO model. The correspond¬ 
ing energy profiles reveal the simultaneous onset of the 
O2-Pd(100) interaction and the SO kinetic energy, and 
therewith underline the decisive role of the SO energy up¬ 
take in ultimately determining the trajectory outcome. 

Quantifying this SO energy uptake in fact reveals that 
considerable amounts of energy are transferred back and 
forth between the O 2 molecule and the SO even on the 
short time scales until dissociation. At the moment when 
the trajectories fulfill the dissociation criterion an aver¬ 
age of 350 meV (considering the sum of kinetic and po¬ 
tential energy of the oscillator) is stored in the SO for 
Tg = 100 K and a low incidence energy of 25meV. For 
the fast molecules with A, = 425 meV this value even 
increases up to 510 meV. In the absence of any further 
dissipation channel, the SO is thus severely overheated. 
As a result, one may expect that the promotion of di¬ 
rect dissociation over “proper” molecular trapping is to 
some degree artificial and that the agreement with the 
experimental sticking curve on the DFT-PBE SO level is 
thereby merely fortuitous. 

This view is indeed supported by the results obtained 
with the DFT-RPBE PES, where, as already indicated, 
dissociation is generally less favored due to the limited 
accessibility of the corresponding transition state. Here, 
the partial damping on account of the surface mobility 
also succeeds in trapping the molecules, but only for a 
slightly extended period of time compared to the ES sit¬ 
uation. This is visible from the extended reflection time 
distribution in Fig. 3f compared to Fig. 3d. The ongo¬ 
ing strong energy exchange between trapped molecules 
and hot SO thus promotes in this case rather a delayed 
reflection than (equally artificially) facilitating dissocia¬ 
tion. The resulting sticking curve at the DFT-RPBE SO 
level shown in Eig. 4 correspondingly still exhibits the 
strong decline with incidence energy. 

Suspecting the lack of further energy dissipation chan¬ 
nels as a major limitation in the present application of the 
SO model we proceed with classical trajectory calcula¬ 
tions that now include a coupling to a bulk thermal bath 
within the GLO approach. At first glance, application of 
the GLO model on the basis of the DFT-PBE PES has an 
almost inconsequential effect on sticking as compared to 
the results obtained within the SO approximation. Eig- 
ure 5 now shows an entirely constant unity sticking co¬ 
efficient independent of both Ei and Tg. Analysis of the 
underlying GLO trajectories nevertheless reveals again 
substantial changes that are fully consistent with the ex¬ 
pectations from the SO analysis. Also at the DET-PBE 
level there is now a contribution of molecular trapping 
to the total sticking, i.e. the GLO succeeds at least to 
some extent in removing the artihcially enhanced direct 
dissociation seen at the SO level. The average amount 
of energy stored in the SO itself at the moment when 
trajectories fulfill the dissociation criterion is now re¬ 
duced to 210 meV (260 meV) at = 25meV (425 meV) 
and Tg = 100 K, as compared to the 350 meV (510 meV) 
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FIG. 5. (Color online) Comparison of the initial sticking prob¬ 
ability So{Ei,Ta) of O 2 on Pd(lOO) as calculated within the 
generalized Langevin oscillator model (GLO) for normal angle 
of incidence and varying incidence energy Ei. Data are shown 
for substrate temperatures of Ts = 100 K (blue circles) and 
400K (red triangles), using either DFT-PBE (dashed lines) 
or DFT-PRBE (dotted lines) energetics. 


found before within the SO model. Owing to the non¬ 
equilibrium nature of O 2 dissociation on Pd(lOO)^^, this 
is still far away from the thermal equilibrium value of a 
surface atom 2 • |A:bTs ~ 25meV even within the limita¬ 
tions of the GLO model, but still suggests more efficient 
surface equilibration. The latter is also reflected in the 
dependence of the trapping contribution to the overall 
sticking on the {Ei,Ts) conditions summarized in Fig. 6. 
As expected for a (largely) equilibrated molecular precur¬ 
sor this contribution depends only weakly on the initial 
incidence energy, but instead sensitively on the substrate 
temperature. As also shown in Fig. 6, the varying ex¬ 
cess kinetic energy of molecules with different Ei is thus 
successfully drained into the heat bath, i.e. the aver¬ 
age energy dissipated into the bath is for every substrate 
temperature precisely by 200 meV higher for molecules 
impinging with Ei = 225 meV as compared to those im¬ 
pinging with Ei = 25 meV. 

On the basis of the DFT-RPBE energetics molecular 
trapping remains the only adsorption mechanism also at 
the GLO level. Allowing for “dissipated” energy to leave 
the SO into the heat bath, however, now suppresses the 
artificial possibility that some fraction of this energy is re¬ 
turned to the adsorbate and thereby induces desorption. 
This removes the delayed reflection times observed within 
the equivalent SO results, cf. Fig. 3f, and yields a largely 
increased sticking coefficient at high incidence energies. 
The resulting DFT-RPBE GLO sticking curves shown in 
Fig. 5 exhibit only a weak decline with Ei. They also ex¬ 
hibit only a weak dependence on substrate temperature. 
This is rather intriguing in view of the predominance of 
the precursor mechanism, which is generally believed to 
fingerprint in form of a sensitive substrate temperature 
dependence. Overall, the DFT-RPBE GLO curve agrees 






Surface temperature, Tg (K) 


FIG. 6. (Color online) Top panel: Temperature dependence of 
the probability of molecular trapping within the GLO model 
and using the DFT-PBE energetics. Bottom panel; Average 
amount of energy dissipated into the GLO heat bath for the 
corresponding trapped trajectories. Data are shown for inci¬ 
dent energies of Ei = 25 (open squares) and 225 meV (closed 
circles). 


therewith now rather well with the experimental data, at 
least similarly well as the DFT-PBE GLO sticking curve, 
cf. Fig. 5. 

This is, in fact, a rather intriguing, if not disturbing re¬ 
sult. At the GLO level both DFT-PBE and DFT-RPBE 
descriptions yield rather similar sticking curves, and this 
despite predicting completely different adsorption mech¬ 
anisms: Predominantly direct dissociation with a small 
contribution of molecular trapping at low substrate tem¬ 
peratures within DFT-PBE versus completely indirect 
dissociation within DFT-RPBE. Compared to the exper¬ 
imental data both functional levels reach indeed the same 
achievements and exhibit the same shortcomings. They 
correctly predict a near-unity sticking over the range of 
studied incidence energies with if at all only a weak de¬ 
pendence on substrate temperature. Yet, both fail to 
reproduce the measured slight increase in sticking with 
Ei, or more precisely they overestimate the experimental 
sticking at low Ei. 

Dissipation seems to affect much more the high inci¬ 
dence energy part of the sticking curves. This view is 
also supported by the robustness of the obtained sticking 
results when varying the parameters entering the GLO 
model as described in Section 11. All in all, these vari¬ 
ations lead only to insignificant variations of So{Ei,Ts) 
within a few percent. In turn, the uncertainties in the 
semi-local DFT description of the actual chemisorption 
well (Eb and E^iss) do critically affect the actual adsorp¬ 
tion mechanism, yet again do not seem to propagate to 
the total sticking curve. In light of the dynamical analy¬ 
ses performed for the FS, SO and GLO models, we sug¬ 
gest that a slightly incorrect description of the entrance 
channel part of the PES can be one reason for the re¬ 
maining small discrepancy - within these effective treat¬ 


ments of dissipation. Presently unaccounted small acti¬ 
vation barriers in this channel for some molecular config¬ 
urations are likely to specifically affect the total sticking 
at low incidence energies and could therewith constitute 
the “missing ingredient” to a fully quantitative agree¬ 
ment with experiment. 

Unless there is a problem in the experimental data, 
the obvious and systematic limitation of these models 
is an equally likely cause: Reducing the phononic fine 
structure of the Pd surface into computationally conve¬ 
nient augmentations of the frozen-surface within the O 2 - 
Pd(lOO) interaction can not, by construction, account for 
the energy exchange with an entire layer of moving sur¬ 
face atoms. Only detailed future studies including all 
these degrees of freedom will allow to completely cap¬ 
ture the influence of surface mobility on the calculated 
sticking curves. 


IV. SUMMARY AND CONCLUSIONS 

We presented a detailed calculation of the initial stick¬ 
ing coefficient of O 2 at Pd(lOO) based on classical trajec¬ 
tory calculations on first-principles 6D potential energy 
surfaces. The specific motivation was to elucidate the 
role of energy dissipation during the exothermic surface 
reaction by augmenting the 6D dynamics with various ef¬ 
fective accounts of surface mobility, as well as to assess in 
how much classical trends in the sticking data reveal the 
underlying adsorption mechanism. To our knowledge for 
the first time, we observe that an account of energy dissi¬ 
pation leads to qualitative changes of the calculated stick¬ 
ing curve as compared to the prevalent frozen-surface ap¬ 
proximation. Gompared to the much more frequently 
studied adsorption of H 2 (also at Pd surfaces^^^^^) this 
might not come as altogether surprising in view of the 
much smaller mass mismatch of oxygen with palladium. 
It is, however, remarkable with respect to the rather sim¬ 
ilar N 2 on W(llO) system® ^ Im)' Intermedi¬ 
ate in its mass ratio compared to H 2 /Pd ~ 

and 02 /Pd ~ exhibits PES character¬ 

istics intriguingly similar to that of O 2 /Pd( 100 ), i.e. 
non-activated paths to dissociation alongside a precur¬ 
sor molecular well. Yet, the N2/W(110) sticking curve 
showed no change when applying effective treatments be¬ 
yond the frozen-surface approach. 

The hierarchical application of FS, SO and GLO mod¬ 
els reveals how the calculated sticking curves sensitively 
respond to details of energy exchange with the substrate. 
Already some degree of energy exchange as provided by 
the SO model enhances the steering of highly accelerated 
impinging O 2 , while “proper” molecular trapping neces¬ 
sitates an account of further bulk dissipation channels as 
in the GLO. The large amount of energy released in the 
exothermic reaction otherwise simply leads to an over¬ 
heating of the SO degree of freedom. Gompared to the 
N 2 /W( 110 ) system, one important difference to note is 
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the considerably deeper molecular well in the correspond¬ 
ing FS-PES description (Eb = —O.SQeV^^ and —0.85eV 
respectively within DFT-RPBE) which leads to a signifi¬ 
cantly increased acceleration of impinging molecules and 
requires more than twice as much energy to be dissipated 
upon reaching the molecular state. The latter may also 
hint towards a different response of the real phonons - 
potentially mitigated by the different phononic proper¬ 
ties of the bare surfaces alone. It might thus comprise a 
key discriminating characteristic, which could be further 
elucidated by more conceptually and computationally de¬ 
manding studies including a more realistic phonon heat 
bath®. 

For two different PES representations based on the 
DFT-PBE and DFT-RPBE functionals we correspond¬ 
ingly obtain at the GLO level a rather satisfying agree¬ 
ment with experimental sticking data^^. This in princi¬ 
ple appealing robustness of the simulation results with 
respect to the uncertainties of semi-local DFT energetics 
does not extend to the level of the underlying adsorption 
mechanism though. The more attractive DFT-PBE ener¬ 
getics predicts a predominantly direct dissociation mech¬ 
anism with some amount of molecular trapping at low 
substrate temperatures. The less attractive DFT-PRBE 
energetics instead predicts adsorption almost exclusively 
via the molecular precursor state. 

Completely different adsorption mechanisms therewith 
lead to rather similar sticking curves that agree equally 
well with the experimental data. Independent of the 
small quantitative discrepancies that remain in either 
case with respect to experiment, this clearly demon¬ 
strates that an unambiguous deduction of the adsorption 
mechanism from the sticking data alone is not feasible 
for this system. Such fingerprinting may work for sim¬ 
pler adsorption systems with smoother potential energy 
surfaces. At the latest for reactions with high degree of 
exothermicity dedicated calculations explicitly account¬ 
ing for high-dimensional potential energy surfaces includ¬ 
ing substrate mobility at best in the form of moving sur¬ 
face atoms are required to establish the mechanistic de¬ 
tails. These details can still be very important not only 
in light of fundamental understanding but also for coarse¬ 
grained micro-kinetic models^^^ . 

Completely different adsorption mechanisms therewith 
lead to rather similar sticking curves that agree equally 
well with the experimental data. Independent of the 
small quantitative discrepancies that remain in either 
case with respect to experiment, this clearly demon¬ 
strates that an unambiguous deduction of the adsorption 
mechanism from the initial sticking data alone is not fea¬ 
sible for this system. Such fingerprinting may work for 
simpler adsorption systems with smoother potential en¬ 
ergy surfaces. At the latest for reactions with high de¬ 
gree of exothermicity dedicated calculations explicitly ac¬ 
counting for high-dimensional potential energy surfaces 
including substrate mobility at best in the form of mov¬ 
ing surface atoms are required to establish the mechanis¬ 
tic details. These details can still be very important not 


only in light of fundamental understanding but also for 
coarse-grained micro-kinetic models^^’^^. 
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